Anomalous Diffusion In Microrheology: A Comparative Study 
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We present a comparative study on two theoretical descriptions of microrheological experiments. 
Using a generalized Langevin equation (GLE), we analyze the origin of the power-law behavior of 
the main properties of a viscoelastic medium. Then, we discuss the equivalence of the GLE with 
a generalized Fokker-Planck equation (GFPE), and how more general GFPE's can be derived from 
a thermo-kinetic formalism. These complementary theories lead to a justification for the physical 
nature of the Hurst exponent of fractional kinetics. Theory is compared with experiments. 
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I. INTRODUCTION 

In the last years, microrheology has become one of the 
most important experimental techniques in soft con- 
densed matter [l|, @, H, 0|- It is a powerful technique 
to determine the viscoelastic properties of complex flu- 
ids at time and length scales complementary to those of 
classical rheological methods [2|, |3j ■ In these complex flu- 
ids, ranging from polymer solutions and colloidal suspen- 
sions to the intracellular medium, the presence of elastic 
forces and molecular motors make anomalous the dynam- 
ics of the Brownian particles 0, [H, @. Hence, anoma- 
lous diffusion becomes a central question that must be 
carefully analyzed due to the presence of confinement 
and finite-size effects related with particle dimensions. 
The dynamics of these particles can be described through 
both, non-Markovian Langevin and Fokker-Planck equa- 
tions @, H, HJ OH [HI- Here, we use and compare two 
different formalisms of anomalous diffusion that in the 
linear force case are equivalent 

The article is organized as follows. In Sec. II we analyze 
the generalized Langevin equation and its equivalence 
with a non-Markovian Fokker-Planck equation. In Sec. 
Ill we introduce the thermokinetic description, show its 
equivalence with the GLE formalism in the linear force 
case, and how it can be generalized to the case of non- 
linear forces. Sec. IV is devoted to the conclusions. 



II. THE LANGEVIN APPROACH TO 
MICRORHEOLOGY 



tion. In these conditions, the motion of the particle could 
be restricted to a small volume and, consequently, its dy- 
namics may manifest confinement and finite-size effects 
modifying the anomalous behavior of the mean square 
displacement (MSD). This last quantity is very impor- 
tant because the viscoelastic properties of the medium 
can be inferred from it through the Stokes-Einstein rela- 
tion 0]. 

In a first approximation, the heat bath can be assimi- 
lated as an effective medium that interacts with the par- 
ticle by means of elastic forces. In first approximation, 
these forces may be represented through a harmonic force 
F/j = — o;, 2 n x, in which uj m is a characteristic frequency 
[4[. Then the dynamics of the particle can be described 
by means of the generalized Langevin equation 



du(t) 
dt 



(3(t-T)u(T)dT + F(t), (1) 



where F(t) is a random force per unit mass and (3{t) is 
a memory function called the friction kernel Q- In the 
Markovian case it takes the form: (3(t) = 0o5(t — r) with 
Po = Qitrja/m the Stokes friction coefficient per mass unit, 
r\ is the viscosity of the solvent and 6(t— r) the Dirac delta 
function. Hence, in this case Eq. fl) recovers its usual 
phenomenological form Q. 

To describe the diffusion of the particle at sufficiently long 
times, t two approximations can be followed. 

a) The Overdamped Case. In this case it is assumed that 
the acceleration term of Eq. ([T]) can be neglected. As a 
consequence the GLE takes the approximate form 



Consider a test Brownian particle of mass m and radius a 
with position x(t) and velocity u(t) = dx(t)/dt. The par- 
ticle performs its Brownian motion through a viscoelastic 
medium which can be made up by a solution of polymers 
or a suspension of particles at sufficiently high concentra- 



0(t - t)u(t)<It + F(t). 



(2) 



It is convenient to stress that in Eq. J2J), the memory 
term involves the velocity of the particle. 
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b) Adiabatic Elimination. In this second approximation, 
one must first solve Eq. © in order to obtain 



u(t) 



X (t-T)x(T)dT+F*(t), 



(3) 



where we have used Laplace transforms, assumed the ini- 
tial velocity of the particle as equal to zero for simplic- 
ity, and defined the memory function x{t) a s the inverse 

Laplace transform x(t) = C - 1 {[s + fa)]' 1 } . This def- 

inition implies that x(^Li s the relaxation function of the 
velocity of the particle [7[. Moreover, we have defined the 
time-scaled random force F*(t) by 



p oo 

F*(t)= / X (t-r)F(t)d7 
Jo 



(4) 



Using now the identity u(t) — dx.(t)/dt = x, one then 
obtains the following generalized Langevin equation for 
the position vector of the particle 

y* oo 

i(t) = -wf n X (t-rMr)dr + F*(t). (5) 
Jo 

In the diffusion regime, t > Pq 1 , Eqs. © and © can 
be considered as the GLE's for the position of the par- 
ticle and then use them to describe anomalous diffusion 
of a Brownian particle in a viscoelastic heat bath. They 
constitute two different models because Eq. © differs 
from © in the fact that it can be obtained through an 
adiabatic elimination of variables, whereas the approxi- 
mated Eq. © has been obtained by neglecting the iner- 
tial term. Physically, the difference lies in the fact that 
the relaxation function of the particle determining the 
time dependence of the MSD, is different. 
Quantitatively, these equations are not simply related 
one to each other, since the memory kernel and the ran- 
dom force are not equivalent due to the scaling intro- 
duced through x(~t)- The fluctuation-dissipation theorem 
(FDT) associated with Eq. © is 



(F(t)F(0)) = 0(t), 



(6) 



where the bracket ( } means the average over noise real- 
izations. The corresponding FDT for the random force 
in Eq. © is given by the expression (F*(t)F*(t')) = 

Jl dr ll dr'x(t - r)*(f - r')/3(r - r'), i- However, in 
the case when the friction kernel is proportional to a Dirac 
delta function, Eq. © defines a thermal (white) noise. 
After performing the adiabatic elimination of variables, 
one obtains the relation (F*(t)F*(0)) ~ exp(-f3 a t) for 
the time-scaled random force. This implies that the posi- 
tion x(i) of the particles satisfies a non-Markovian equa- 
tion with a exponentially decaying time correlation of the 



random force, Q. Only in the Markovian case Eqs. © 
and (J5J) are equal. 

At certain time scales, anomalous diffusion is character- 
ized by a power-law behavior of the MSD of the particle 
as a function of time. Within a Langevin description, this 
behavior is a consequence of the statistical properties of 
the random force, that is, of the heat bath as follows 
from the relations already obtained. In particular, by us- 
ing Eqs. ([2]) and ©, it can be shown that anomalous 
diffusion is a consequence of a fractionary Gaussian noise 
(FGN) that satisfies the FDT [l(J 



(F(t)F(0)) = (k B T/m)f3 2 2H(2H - l)\t\ 



2H-2 



(7) 



where the exponent must satisfy the relation: 1/2 < H < 
1 , the factor /3q gives the correct dimensions and we have 
introduced the dimensionless time t — t ~ 1 t, with to a 
characteristic time. In particular, to can be taken as the 
characteristic diffusion time tjj = a 2 !?^ 1 , with Do = 
fcsT/(m/3o). Hence, using Eqs. © and ([7]) one may 
calculate the time dependence of the MSD by using the 
method of Laplace transforms. Here, it is convenient to 
write ([2j in the dimensionless form 



x(t) = -a j3(t- f)u(f)df + F(i). 



(8) 



where (3(t) = \t\ 2H ~ 2 and F = oT^lq^F are dimension- 
less quantities in accordance with Eqs. (J7J and ©, and 
a = 2H(2H - l)P 2 Lu^ t 2 . In writing © we have used t 
and the dimensionless position vector x = a _1 x. 
Now, by Laplace transforming Eq. ([8]), one may obtain 
the expression for the transformed velocity u(s), which in 
turn can be substituted into the definition sx(s) = u(s) 
with x(0) = for simplicity. After solving for x(s) one 
obtains 



x( S ) 



- 1 T~ 1 (2H - 1) 



,2ff-2 



1 



l T- 1 (2H-l)s 



2H-2 



(9) 



where T(2H — 2)s 1 2H and F(s) are the Laplace trans- 
forms of (3 and F, respectively. The inverse Laplace 
transform of Eq. © is given in terms of the convo- 
lution of the noise F and the Mittag-Leffler polyno- 
mial zl~ 2H E 2 -2H,2-2H (—z 2 ~ 2H ), with the scaled vari- 
able z* = [aT{2H - i)]" 1 /(2-2H) 5; ^ 
The expression for the MSD can now be obtained by cal- 
culating (x 2 ){t) in terms of the inverse Laplace trans- 
form of Eq. ©, using the FDT © and performing the 
corresponding integrals. The general expression involv- 
ing Mittag-Leffler polynomials is complicated, however at 
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times t satisfying (3 <C t ~ to, one obtains the following 
short-time power law behavior 



(x^w-^r^K) (t/to) 



2-2H 



(10) 



where some constants have been dropped for the sake of 
simplicity. 

In Figure [TJ we compare the MSD given by Eq. 
(flT)]) (dash-dotted lines) with experiments (symbols) of 
0.965/im diameter latex microspheres imbedded in F- 
actin solutions of increasing concentrations, [3j. As ex- 
pected, the agreement between theory and experiments is 
good for short times Pq 1 <C t ~ to with to ~ td playing 
the role of a crossover time. If we take f3o ~ 10 6 s _1 then 
td ~ 10~ 2 s in typical situations. 

The MSD can be used to obtain the microrheological 
properties of the medium by means of the generalized 
Stokes-Einstein relation 0. From Eq. © in the zero 
inertia approximation, one may show that the complex 
shear modulus of the medium G"(cu) satisfies the rela- 
tion 



G"{uj) ~ (k B T/muj(x 2 )(uj)) , 



(11) 



where u) is the frequency. Now, by taking the Fourier 
transform of Eq. (flU|) and substituting the result into 
(fTTj) . one obtains the scaling relation G"{tu) ~ (touj) 2 ^ 2H ■ 
From this analysis it follows that the power-law behavior 
of the MSD of the particles and the related complex shear 
modulus of the effective medium, are a consequence of the 
fractionary Gaussian noise of the GLE Notice that 
Eq. ([5]) can also be used to describe anomalous diffusion, 
however this will be explained after Eq. (fl~5]) . below. 
The Generalized Fokker-Planck Equation. An equivalent 
approach to study anomalous diffusion in a viscoelastic 
bath can be performed by means of a generalized Fokker- 
Planck equation . 

By using the Laplace transform technique, in Ref. 0| it 
has been shown that the GLE ([1]) is equivalent to the fol- 
lowing non-Markovian Fokker-Planck equation contain- 
ing time dependent coefficients 



I 8£ 
dt J dx 



(u/) 



w 2 (f)x 



ou ou 



dx 



d ~ 
Ou 



u/ + 



k B Tdf 



du 



(12) 



Here, /(x, xo; u, Uo, t) is the probability distribution 
function depending on the instantaneous position r 
and velocity u of the Brownian particle. ro and 
Uo are the corresponding initial conditions. In this 
equation, memory effects are incorporated through the 
time dependent coefficients /3(t) and ui(t) and ip(t) — 
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FIG. 1: MSD vrs time. Comparison between theory (lines) 
and experiments with actin filament networks (symbols) 
(taken from Ref. Q|). The dashed line corresponds to the 
power law behavior valid for long times, Eq. (|16p with 
to = 2 • 10~ 3 s, u m ~ 2.6 ■ 10 2 s"\ and b ~ 1 ■ 10" 2 . The 
dotted lines correspond to a stretched exponential behavior 
Eq. ((2TJ>, with lu^/uj^ ~ 5.61 ■ 1CT 4 , 2Lj^B 2 t /apo ~ 35 and 
2/So/Qto^m =3/4. The solid lines correspond to the nonlinear 
model Eq. (|28[1 . The values of the parameters are indicated 
in the text. The dash-dotted line corresponds to Eq. (|10[) 
with H = 5/8. The numbers 24, 63 and 124 correspond to 
increasing concentrations of actin filaments, 



{ksT /muj 2 n )\uJ 2 {t) — (J^. They are given in terms 
of combinations of the relaxation functions Xu{t) — 
m(3/c B T)- 1 (x(i) -u ) and X *(i) = (xg)" 1 (x(i) • x ), and 
their explicit form is given in Ref. [7|. It is worth to 
stress here that an equation similar to (fT2|) has also been 
obtained by using projector operator techniques in Q. 
The result expressed in Eq. (|12p is important, because 
it implies that generalized Fokker-Planck equations in- 
corporating memory effects through memory functions 
constitute models not simply related with Eq. (fTJ), [l|, 0]. 
Passing from one model to the other implies an approxi- 
mation that imposes conditions on the rate of change of 
the fields involved in the description [1, [ll| . 
At times t Po 1 , from Eq. (fl"2")) one may derive a 
generalized Smoluchowski equation for the mass density 
p(x, t) = m J fdu. To derive it, one may follow the adi- 
abatic elimination of variables procedure by calculating 
the evolution equations for the first three moments of 
/(x, u,i). Then, after imposing the condition of large 
times t 3> Pq 1 i a reduced description is obtained in terms 
of an evolution equation for p. 

Averaging Eq. (fl"2")) over u assuming that the currents 
vanish at the boundaries, we obtain the continuity equa- 
tion 



d P /dt = -v • (pv), 



(13) 
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where v(x, t) — mp^ 1 J ufdu is the average velocity field. 
At long times one obtains the following constitutive rela- 
tion for the diffusion current lUSl 



-p~ 1 u 2 x.p-DVp, 



(14) 



where we have defined the effective diffusion coefficient 
D(t) — (k B T / muJ^) ft -1 Co 2 and used the definition of 
ip(t). Substitution of Eq. (QJJ into §T3§ leads to 

dp/dt = D(t) [V 2 p - (mu 2 Jk B T) V • (px)] , (15) 

which constitutes the generalized Smoluchowski equation 
(GSE) for p. An equation similar to fT5]) can also be ob- 
tained from Eq. ([T]) by using Laplace transforms. In 
Ref. the authors found that (m/k B T)D(t) is pro- 

portional to — dln\R(t)\ /dt, with R(t) = x(t)/x(to) an d 
x(t) — (x ■ xo)(i) the relaxation function of the position 
of the particle. If one considers the case of FGN, the 
relaxation function R(t) is given in terms of a Mittag- 
LcfHcr polynomial leading to a complicated form for the 
effective diffusion coefficient D{t), By evaluating 

numerically the obtained relation, it can be shown that 
at short times the function follows a power law behavior 
with an exponent equal to that calculated by simply tak- 
ing the time derivative of Eq. (JTUJ). At larger times, a 
crossover arises in which the relaxation function can be 



assumed to decay as the power law R(t) 



Using 



the definition (x 2 ) — J x 2 pdx, Eq. (|15p and assuming 
/3 _1 u) 2 = /3 _1 cij 2 n 7(t), with 7 dimensionless, one may cal- 
culate the following expression for the MSD 



(x 2 ) ~ 3(k B T/muj 2 m ) 1 - (i/i„) 



,-2b 



(16) 



At times t ^> t , Eq. (fT())) can be used to fit the experi- 
mental results, as shown by the dashed line in Figure [TJ 
See the discussion after Eq. (fit))) . 

This analysis shows that, in the linear force case the GLE 
with FGN offers a good explanation of the subdiffusion 
observed in microrheological experiments. A different 
approach in order to have deeper understanding of the 
physical mechanisms responsible for subdiffusion in a vis- 
coelastic fluid will be presented in the following section. 



THERMOKINETIC-DESCRIPTION OF 
ANOMALOUS DIFFUSION 



III. 



The Fokker-Planck type Eqs. (|T2j) and (|T5j) for the prob- 
ability distribution function, admit a different interpreta- 
tion in terms of kinetic equations. This rcintcrpretation 
enriches and simplifies the description of anomalous dif- 
fusion, because allows one to formulate phcnomcnological 



models for the coefficients $ and uj which not do depend 
on the a priori election of the statistical properties of the 
random force of Eq. fT} . 

Fokker-Planck type kinetic equations for different phys- 
ical phenomena can be obtained by using two general 
principles. The first one is the conservation of the prob- 
ability, expressed by the equation 



9/(r,i)/9t = -V r •[/(£, *)V 7 ], 



(17) 



where T represents the set of variables necessary to de- 
scribe the state of the system at a mesoscopic level. Here, 
f(Hi represents a diffusion current in T-space and 
Vr the corresponding gradient operator. The second 
principle is the so-called generalized Gibbs entropy pos- 
tulate 0, Uli 

P As(t) = -k B J f(T,t)ln\f(T,t)/fi eg (T)\dT, (18) 

where As is the difference of the specific entropy with 
respect to a local equilibrium reference state, character- 
ized through the probability distribution fi eq {Q- The 
assumption of local equilibrium allows one to calculate 
/ierj(r) by using the techniques of equilibrium statistical 
mechanics. 

To obtain the evolution equation for /(r, i), we may fol- 
low the rules of mesoscopic nonequilibrium thermody- 
namics (MNET) [I, HI]. Schematically, by taking the 
time derivative of Eq. (fT5|) and inserting (fTT)) . after inte- 
grating by parts assuming the usual boundary conditions, 
one obtains the entropy production o~(t) 



<r® = ~ J /(r,t)V 7 .V r 



k B T 



rn 



m|//./WI 



dT. (19) 



In a similar way as in linear irreversible thermodynam- 
ics [l5j . we can assume linear relationships between the 
forces {ksT /m)h\\f / fi eq {T)\ and the conjugated cur- 
rents /V T . Using them, one finally arrives at the multi- 
variate Fokker-Planck equation 



df/dt - V r • e(t) [(k B T/m)V r f - (/X)] 



(20) 



where the generalized force X(r) is defined in the usual 
form: X = -^TVr ln/; e(? , [9J. Here, memory effects are 
incorporated through the time dependent Onsager coef- 
ficients introduced by the linear coupling between forces 
and fluxes, and contained in £(i), [J, Il3j |. The tenso- 

rial character of £(i) accounts for the anisotropy of the 
medium. 

If we consider a Brownian particle under the action of 
a linear force and whose mesoscopic state is determined 
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by its instantaneous velocity u and position x, Eq. (|2TJ)) 
takes the form (fT2|) . In this case, the local equilibrium dis- 
tribution is fi e q — foexp [— (fcsT/2m) (u 2 + 20)] , where 
</>(x) is the potential associated with the harmonic force, 
fo is a normalization factor and we have assumed that 
in T-space entropic forces are coupled to the currents in- 
dependently from those arising from an energy potential, 
reflecting its different physical origin [l3T |. 
At larger times, we may assume that the state of 
the particle is only determined by its position vec- 
tor, then the local equilibrium distribution is pi eq = 
poexp [— {ksT / m)4>{'x)\, where the potential 0(x) is ar- 
bitrary. By following the procedure indicated above, the 
following GSE can be derived 

dp/dt = D(t) {V 2 p - (m/k B T) V • [pX(x)]} , (21) 

where we have assumed an isotropic medium when writ- 
ing the scalar effective diffusion coefficient D(t) and de- 
fined X = — V0(x). Hence, the thermokinetic formalism 
allows one to derive generalized Fokker-Planck equations 
of the form (|2~lj) which contain nonlinear forces. No as- 
sumptions on the statistical nature of the random force 
have been done until now. This is an important difference 
with respect to the GLE description, because the method 
presented in the second section gives exact Fokker-Planck 
equations only in the case of linear forces. 
During its motion through the viscoelastic medium, the 
Brownian particle induces perturbations on the velocity 
field of the host fluid surrounding it. These perturba- 
tions propagate and are reflected by the local boundaries 
which can be made up by a polymer network or another 
suspended particles. As a consequence, they modify the 
velocity field of the fluid around the particle in a later 
time, introducing memory effects. In this form, hydro- 
dynamic interactions change the distribution of stresses 
used to calculate the force over the surface of the particle. 
In first approximation, this force is characterized by the 
effective mobility coefficient -f(t) = (m/kBT)D(t), in the 
form [J] 



7(f) - P^a^t). 



(22) 



Here, 7(f) is a dimensionless function of time and a — 
1 + B 1 a/y is the mentioned correction 0,[I3,[3- The co- 
efficient Bi depends on the nature of the boundary (solid 
wall, membrane, polymer network or a cage formed by 
surrounding particles) and y represents a characteristic 
length of the medium. 

By using the dimensionless variables previously intro- 
duced and defining the time scaling r(t) = J ^(t)dt, from 




o Experiments 

Theory, Markoviatt approximation 
— Theory; non-Markovian approximation 



FIG. 2: MSD vrs time. Comparison between theory (lines) 
and experiments of chromatin diffusion in living cells. The 
dotted line corresponds to a linear model in the Markovian 
approximation. The data were taken from Ref. [f| (circles) 
and the figure from [20| . 



(j2Tj) one obtains 

dp/dr = a (u^h/Po) V 2 p - a (F i //3 a) V • ( P x) , 

(23) 

where lo\ = fc^T/ma 2 and Fo is the magnitude of the 
force per unit mass X due, for instance, to the polymer 
network. 

The time dependence of r can be obtained by using 
the evolution equation for the time correlation function 
X (t) = (x-x )(f) (see, Eq. [15]), Q 



d X (t)/dt = ~a {F t /f3 a) 7(f) (X • x )(t). 



(24) 



Eq. (|24p involves a complete hierarchy for the moments 
of the distribution p. However, in the linear force case it 
takes a closed form leading to the relation 



7(f) = din 



R -fh/ a t u> m ( 25 ) 

2 x with ux^ = a~ 1 Fr j and 
R(t) = X(t)/x(0)- Using the definition and Eq. 



where we have used X 



one gets r(i) = In 



These relations imply: 



D = D, see the discussion after Eq. (fl"5|) . 
At short times, we can infer the time dependence of 
T(t). Expanding the logarithm in its argument around 
the unity: In \R~ P°/ at °"l | ~ R -f3 / a t Q u 2 m _ 1 + q(R 2 ) 
and taking into account that R(t) must be an even func- 
tion of time, a first order expansion leads to: i? _1 (t) ~ 
1 + B^t 2 + 0(t A ). Thus, we obtain the approximate ex- 
pression 



(26) 



where the parameter B2 oc B^ in general depends on the 
characteristic length of the medium 4]. Important to 
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FIG. 3: MSD vrs time. Comparison between experiments 
(circles) and theory (solid line) for the case of actin fila- 
ment networks. The data were taken from Ref. 0|. The 
theoretical model corresponds to Eq. (128[) with {x 2 )o ~ 



4.4 ■ 10"V"i 2 , Do ~ 9.3 ■ 10" 5 /i"i 2 s"\ /3o ~ 5.5 • 10V 
F ~ 3.8 ■ 10" 

Q ~ 1. 



Urns' , t ~ 0.1s, A ~ 1.5 • 10~>m, B 2 



60, 



notice is that the exponent depends on /3q, to, uj m and 
a. Eq. (j2"rj)) is consistent with a stretched exponential 
relaxation of the correlations [l^|. Hence, in the linear 
case, Eq. Q gives the MSD 



(x 2 )(r) = 3(<4/<«&) 



1 



,-2acj^t O T(t)/0 o 



(27) 



which at short times yields (x 2 ) ~ 6 J B 2 ^? T P /3o/QtoW 
Comparison with Eq. 



to 



(3§uj~ 4 and H = 1 



/3o* 

TU|) leads to the identifications 
Po/atouj 1 ^. This expression for 



constitutes the central result of this work. It shows 
that the so-called Hurst exponent is a consequence of hy- 
drodynamic and elastic interactions between the particle 
and the viscoelastic medium. 

This model explains the behavior of the MSD in the case 
of constrained diffusion Q. Fig. [2] shows a compari- 
son between theory (solid line) and experiments of con- 
strained diffusion of chromatin in living cells [f|. The 
dashed line corresponds to a similar model in the Marko- 
vian case. As expected, the short time behavior of the 
MSD is subdiffusive. These effects are also present in Fig. 
Q] (dotted lines). At short times the behavior of MSD fol- 
lows a power law on the scaled time t and saturates at 
long times, where confinement becomes significant. A 
more detailed model is reported in j4j. 
A nonlinear model can also be discussed in terms of Eq. 
(f2"Tj) . If we assume that the force exerted on the particle 
by the medium is of the form cos [A _1 (x)(t)] , with A 
the average distance between the particle and the local 
boundaries, it can be shown that 



(x 2 )(t) ~ (x 2 ) + 2D r{t) + 4A 2 tanh 



FqDq 

2\k B T 



r(t) 



where (x 2 )o is some initial value. The behavior of the 
MSD (j28|) as a function of time (solid line) is shown in 
Figs. 1 and 3, and compared with experiments (symbols). 
The plateau is a signature of the existence of cage effects, 
which in our model are related with the maximum value 
of the elastic force X. The agreement between theory 
and experiments is good. 



the values of 
2 Dn ~ 1.3 • 



For Fig. 1 in the case with label 164 
the parameters are: {x 2 )o ~ 5 • 10 -7 /ito 
1CtV« 2 s~\ A) ~ 5.5 • 10 9 s-\ F ~ 1 ■ JLU "(ims 
t ~ 0.1s, A ~ 1.1 • l(T 2 ^m, B 2 ~ 7.2, a ~ 1 and 
we have defined w 2 ^ ~ FqX" 1 . Finally, it is interesting 
to notice that the effective mobility 7(4) may be inter- 
preted phenomenologically as the inverse Laplace trans- 
form: Cr 1 [uj/G"{oj)]. Taking into account the relation 
between the MSD and the compliance J: (x 2 )(t) ~ J(t), 
[3], one obtains j(t) = /T 1 [uj 2 J{uj)]. 



(28) 



IV. CONCLUSIONS 



By following two different approaches, in this article we 
have described the anomalous diffusion of a Brownian 
particle in microrheological experiments. In the linear 
force case, we have shown that the generalized Langevin 
description with fractionary Gaussian noise explains well 
the behavior of the microrheological properties of the vis- 
coelastic medium. We have also shown the equivalence of 
this description with generalized Fokker-Planck equations 
having time dependent coefficients. The comparison with 
a thermokinetic formalism gives a plausible justification 
for the physical nature of the exponent H = l—fto/atoU)^ 
used in the GLE description. The thermokinetic formal- 
ism constitutes a powerful generalization of the theory 
allowing the formulation of non-Markovian models even 
in the case of non- linear forces. The agreement between 
theory and experiment is good. 
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